Rate effects on layering of a confined linear alkane 
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We perform drainage experiments of a linear alkane fluid (n-hexadecane) down to molecular 
thicknesses, and focus on the role played by the confinement rate. We show that molecular layering 
is strongly influenced by the velocity at which the confining walls are approached: under high 
enough shear rates, the confined medium behaves as a structureless liquid of enhanced viscosity for 
film thickness below ~10 nm. Our results also lead us to conclude that a rapidly confined film can 
be quenched in a met ast able disordered state, which might be related with recent intriguing results 
on the shear properties of confined films produced at different rates [Zhu and Granick, Phys. Rev. 
Lett. 93, 096101 (2004)]. 

PACS numbers: 81.40.Pq, 68.35.Af, 83.50.-v 



When liquids are confined to molecular thicknesses, 
their properties deviate markedly from those of the bulk 
state [l|, l2[ . If two solid surfaces are approached at dis- 
tances of a few nanometers in a liquid, so-called struc- 
tural forces are observed, which are associated with the 
ordering of molecules into layers parallel to the confining 
walls 11, 3, 4, Understanding how such a struc- 

turation affects the shear response of confined liquids has 
been the scope of a large number of studies, both numer- 
ical [5|, |y, [Z| and experimental. The latter have adressed 
the issue of flow in nanometer-thick films usin g th e Sur- 
face Force Apparatus (SFA)0^,[li,[l3,[il,M[3 or 
Atomic Force Microscopy [18|, LL9| . 

Most of the works mentioned above focus on the lu- 
brication or viscoelastic properties of equilibrated films, 
i.e. which have reached their equilibrium layered struc- 
ture before shear begins. However, in the course of a 
recent debate on reassessement of SFA measurements 
[23I . I24I I25I ]. it has been suggested that the velocity at 
which a liquid is brought to molecular thickness may 
strongly influence its shear response [23I . [24|. In their 



work, Zhu and Granick 23|, l24| rind that thin films 



formed by rapid confinement exhibit a much higher ef- 
fective viscosity than those produced quasistatically. The 
authors attribute this effect to the fact that rapid con- 
finement yields less structured films. 

This raises two questions: (i) how does molecular lay- 
ering depend on the confinement rate, and (ii) if a more 
disordered film is produced by rapid confinement, what 
is the stability of such a "mechanically quenched" state 
? 

These questions have motivated us to revisit drainage 
experiments, as performed initially by Chan and Horn 
[ill , and to study how the thinning of n-hexadecane con- 
fined between mica surfaces is affected by the velocity 
at which the surfaces are approached. We show here 
that layering is indeed most sensitive to the confine- 
ment velocity: at high enough shear rates, structuring 
is completely hindered. For thicknesses in the range 3- 
10 nm the drainage dynamics is akin to that of a liquid 



whose effective viscosity increases with the level of con- 
finement. Furthermore, our results lead us to conclude 
that a non-structured film of hexadecane, obtained by 
quenching rapidly down to ~ 2 nm, is metastable and 
relaxes towards the layered equilibrium configuration via 
a nucleation/growth process. 

The experiments were performed using a recently de- 
veloped surface force apparatu s |26f - As in previous 
versions of the instrument [13, l28l. I29I]. two atomically 
smooth mica sheets are mounted in a crossed-cylinder ge- 
ometry (radius of curvature R ~ 1 cm), and the surfaces 
are immersed in the liquid under study. Our appara- 
tus allows for independent measurements of the normal 
force F (by means of a capacitive load cell of stiffness 
31000 N.m -1 ) and intersurface distance d (using multi- 
ple beam interferometry) at a rate on the order of 30 Hz. 
This enables us to obtain F(d) profiles during drainage 
of the fluid, while the mica surfaces are approached by 
driving the remote point of the loading spring at a pre- 
scribed velocity V in the range 0.05-20 nm.s -1 . The mica 
sheets were prepared according to the following proto- 
cole. Muscovite mica plates (JBG-Metafix, France) were 
cleaved down to a thickness of ~ 10 /im, cut into 1 cm 2 
samples by means of surgical scissors, and coated on one 
side with a 40nm-thick thermally evaporated silver layer. 
The sheets were fixed, silver side down, onto cylidrical 
glass lenses, using a UV setting glue (NOA 81, Norland. 
UV curing was followed by thermal aging at 50° C for 
12 hours in order to suppress visco-elasticity of the glue 
layer). The mica sheets were placed on the lenses so that 
their cristallographic axis be aligned. Prior to the exper- 
iments, each mica sample was recleaved using adhesive 
tape [30], and mounted in the apparatus so that the re- 
gion of closest distance between the lenses be free of steps. 
The surfaces were brought into contact under an argon 
atmosphere, and the total mica thickness was deduced 
from the position of the fringes of equal chromatic or- 
der (FECO) using the multilayer matrix method [26j.[3lj. 
The surfaces were then separated and the thickness of 
each mica sheet determined using the same method. A 



2 



drop of liquid (~30 fjL) was finally injected between 
the surfaces, a beaker containing P2O5 was placed in- 
side the apparatus which was then sealed and left for 
thermal equilibration for about 12 hours before measure- 
ments began. The liquid used was the linear alkane n- 
hexadecane (Sigma- Aldrich, >99%). The product was 
filtered through a 0.2 /am membrane immediately before 
injection. All the experiments reported below have been 
performed at T = 24 ± 0.02°C. 

We present results obtained in the following conditions. 
Starting from an initial distance of 20-30 nm, the remote 
point of the normal loading spring was driven at a given 
speed V, until a normal force on the order of 500 /iN was 
reached. Elastic flattening of the surfaces was observed 
for loads larger than 100 /iN, and F =500 /iN yields a 
normal pressure of about 500 kPa. 

Figure [T] shows three F(d) profiles measured at V = 
0.05, 1, 5 and 10 nm.s -1 . Layering of hexadecane is 
clearly observed at V = 0.05 nm.s -1 : for distances d < 3 
nm, the film thickness decreases by steps of 4-5 A as the 
force increases, in agreement with previous studies [!,[l(j. 
Such a steplike thinning of the fluid totally disappears 
for V > 1 nm.s -1 to yield a smooth monotonic repulsive 
profile (see Fig. [I]). Besides, we note that for a given 
film thickness, the normal force is all the larger as the 
confining velocity is high. 

So, for high enough approach velocities, layering is 
dynamically hindered, and the confined medium stays 
structureless. 
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FIG. 1: Force vs distance curves measured during approach 
of the mica surfaces at V = 0.05 nm.s -1 (full line); V — 1 
nm.s -1 (full line with o markers); V = 5 nm.s -1 (full line with 
• markers); and V = 10 nm.s -1 (full line with A markers). 
Inset: close-up view of the same data set. Arrows indicate 
the first layering transition detected, between d ~3 and 2.5 
nm. 

Now, at large V, hydrodynamic forces may account for 
a non- negligible part of the measured repulsive forces. 
On fig. [2j we compare the force-distance curve measured 



at V = 10 nm.s -1 with that corresponding to the squeeze 
flow of hexadecane, using the bulk viscosity 77 = 3.5 x 
10 -3 Pa.s. The hydrodynamic force, assuming no-slip 
boundary conditions, is given by [111 ]: 

F H = ^R 2 v^ (1) 

where R = 1 cm is the radius of curvature of the cylin- 
drical lenses, and d is the actual approach velocity of the 
surfaces. Comparison of Fn(d) and F(d) is performed 
for F < 100 /iN, to ensure that the radius R is not af- 
fected by flattening of the surfaces. Fig. [2] shows that for 
d >10 nm, the repulsive force can be ascribed to hydro- 
dynamic flow of the bulk liquid, whereas for d < 10 nm, 
the measured forces are larger than Fh- A modification 
of eq.(pQ), as suggested by Chan and Horn [ll|, to account 
for shifted no-slip boundary conditions — i.e. replacing 
d by d — 2d si where d s is the thickness of an immobile 
layer near each wall — does not allow either to fit the 
experimental F(d) profiles at short distances. 




-1 10' 5 I i 1 1 1 1 

5 10 15 20 25 30 

d (nm) 



FIG. 2: Force vs distance curve measured during approach 
of the mica surfaces at V = 10 nm.s -1 (full line). Hydrody- 
namic force Fh calculated using expression (pQ) and a no-slip 
boundary condition (full line with •). Fh calculated using a 
no-slip plane shifted inside the gap by d s =0.8 nm (full line 
with o). 

This points to the existence of a thickness regime in 
which the non-structured confined fluid does not ex- 
hibit bulklike flow properties. Tentatively, we use ex- 
pression (pQ) to estimate an effective viscosity of confined 
hexadecane from the experimental F(d) data: 7/ e ff — 
Fd/(67rR 2 d). This is plotted on Fig. [3] as a function of 
thickness, for three driving velocities. It is seen that com- 
parable profiles for rj e fi(d) are obtained, independently of 
V, which gives support to our analysis in terms of an ef- 
fective viscosity. Moreover, it appears that ^ e ff starts to 
deviate noticeably from the bulk value at short distances, 
and increases by more than one order of magnitude as d 
goes roughly from 10 to 3 nm. 
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This is in apparent contradiction with the results of 
Chan and Horn [ll|, who found that, assuming the exis- 
tence of two immobile layers on each wall, bulk viscosity 
could account for the fluid drainage down to thicknesses 
of 2-3 nm. However, it is important to note here that 
our apparatus has a stiffness which is more than two 
orders of magnitude larger than that used in reference 
From the force balance Fh = K(Vt — [d — dt=o])> 
where the right-hand-side term is the restoring force of 
the cantilever spring of stiffness K, it follows that for a 
given distance d and driving velocity V, the lower 
the lower the approach velocity <i, as well as the max- 



imum shear rate to which the fluid is submitted [111 ] : 
7max ~ yR/d x d/d (7 reaches its maximum value at a 
lateral distance of a few microns away from the point of 
closest approch). In the inset of Fig 131 we have plotted 
7max(^) for V = 10 nm.s -1 . It appears that r] e Q becomes 
larger than ?7buik when 7 max > 10 3 s -1 . An analysis of the 
data of Chan and Horn in order to evaluate 7 ma x shows 
that, in their experiments, 7 ma x < 10 3 s -1 , and that its 
maximum value is reached at distances d ~ 3 nm (see 
inset of Fig(3]). Their study and ours are actually com- 
patible, and the present work highlights the advantage 
of a high stiffness SFA, which gives access to a regime 
which could not be investigated in previous studies. 

Furthermore, from the shear rate value of 10 3 s -1 at 
which rjefi > r^buik, we estimate the relaxation time of 
hexadecane molecules to be l/7max — 10 -3 s when d < 
10 nm, i.e. many orders of magnitude larger than the 
Rouse time of bulk Ci6 chains [32|. This is in qualitative 
agreement with previous studies on relaxation in confined 
fluids Q. 



Coming back to Fig. [TJ it is clear that under a normal 
force of 500 /iN, the thickness of the film is 1.65 nm when 
V = 0.05 nm.s -1 and 1.9 nm when V = 10 nm.s -1 . In or- 
der to determine whether the thicker, rapidly quenched, 
film evolves with time, we have performed the following 
experiment: at V = 10 nm.s -1 , hexadecane was confined 
until a normal force of 500 /iN was reached. The load 
was then maintained at this value by means of a feed- 
back control loop, and we monitored the time evolution 
of the film thickness d(t), which is plotted on Fig. [4j The 
d(t) curve exhibits three important features: (i) once the 
force setpoint has been reached (t = on Fig. [4j), the 
film thickness stays constant at 1.9 nm for ~20 s; (ii) for 
t > 20 s, d decreases linearly with time (d = 3 pm.s -1 ), 
until it reaches a constant value of 1.65 nm; (iii) the final 
value of d corresponds to the thickness of 4 molecular 
layers obtained under quasi-static loading. 
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FIG. 3: Effective viscosity as a function of film thickness for 
V — 2 nm.s -1 (line with • markers), V — 5 nm.s -1 (full line), 
and V = 10 nm.s -1 (full line with o markers). The dashed line 
indicates the bulk value ?7 = 3.5xl0 -3 Pa.s. Inset: maximum 
shear rate jmax vs distance for V = 10 nm.s -1 (full line) in 
our experiment, and that deduced from Fig. 9 of Ref. [Ill ], for 
which V = 16 nm.s -1 (line with o markers). 



FIG. 4: Thickness vs time for two films confined at driving 
velocity V = 10 nm.s -1 (black line), and V = 1 nm.s -1 (red 
line). Time t = corresponds to the point were loading is 
stopped and a normal force of 500 /iN is maintained, t < 
corresponds to the end of the loading phase. The thick lines 
are guides for the eye. 

We propose the following picture to interpret such a 
behaviour. The abrupt stop of the thickness decrease, 
together with its quasi-linear shape, lead us to rule out 
relaxation via diffusive- like evacuation of free volume to- 
wards the contact edges, which would yield a gradual self- 
deceleration. On the contrary, such a behavior appears 
consistent with a nucleation/growth process by which the 
rapidly confined film performs a transition from a higher 
volume, metastable disordered state, to a lower volume 
stable layered structure. A similar mechanism has al- 
ready been evidenced in the layer- by- layer collapse of 
confined liquids [ljj. In this context, one is tempted to 
identify the initial 20 s "incubation" time with a nucle- 
ation delay. However, since our signal is an average over 
a (5 /im) 2 zone close to the contact center, we cannot 
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exclude that it might correspond to the time required for 
the growing layered patch to reach a large enough radius 
(comparable to the mica thickness ~ 3/im) for the local 
bending of the mica sheets to be detectable. 

The above picture is further supported by the follow- 
ing control experiment: hexadecane was confined at a 
much lower rate, namely V = 1 nm.s -1 , and d(t) moni- 
tored under F =500 /xN. It is seen on Fig. [4] that, once 
the force setpoint is reached, the film thickness decreases 
from d = 1.75 nm at the same velocity and until the 
same final value as those observed in the higher rate ex- 
periment. This is consistent with the growth of a stable 
phase which, under a given applied pressure, occurs at a 
velocity which is independent of the rate at which the film 
has been produced. Moreover, we note on Fig. [4] that no 
"nucleation phase" is observed for V = 1 nm.s -1 , and 
that d decreases with time during the loading phase, at 
a velocity which is about twice that observed under con- 
stant load. This suggests that nucleation of the stable 
phase has occured before reaching the load setpoint, and 
that propagation, helped by the slowly increasing pres- 
sure, takes place at a larger velocity during loading. 

In conclusion, we have performed drainage experi- 
ments which unambiguously show that layering of a con- 
fined liquid is hindered when the medium is submitted 
to high enough shear rates. Taking advantage of a newly 
developed Surface Force Apparatus, we evidence that 
the effective viscosity of hexadecane may deviate from 
its bulk value at thicknesses on the order of 10 nm in 
situations of fast out-of-equilibrium confinement. More- 
over, our results indicate that a disordered film formed 
by a rapid mechanical quench is in a met ast able state. 
Whether the formation of such metastable states is at 
the origin of the intriguing results of Zhu and Granick 
[23L 0] is still to be clarified. This calls for further ex- 
periments investigating how the yield stress build-up and 
the steady state shear response of molecularly thin films 
are affected by the confinement rate. 

The author wish to thank C. Caroli and T. Baumberger 
for fruitful discussions. 
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